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Abstract. We present the results of ISO SWS and LWS grating 
scans towards the three brightest members of the BD+40°4124 
group in the infrared: BD+40°4124 (B2Ve), LkHa 224 (A7e) 
and the embedded source LkHa 225. Emission from the pure 
rotational lines of H2, from ro- vibrational transitions of CO, 
from PAHs, from H I recombination lines and from the infrared 
fine structure hues of [Fe ll], [Si ll], [S l], [O l], [O III] and [C ll] 
was detected. These emission lines arise in the combination of 
a low-density (« 10^ cm^'^) H II region with a clumpy PDR in 
the case of BDh-40°4124. The lower transitions of the infrared 
Hi lines observed in BDh-40°4124 are optically thick; most 
likely they arise in either a dense wind or a circumstellar disk. 
This same region is also responsible for the optical H I lines 
and the radio continuum emission. In the lines of sight towards 
LkHa 224 and LkHa 225, the observed emission lines arise in 
a non-dissociative shock produced by a slow (« 20 km s^^) 
outflow arising from LkHa 225. Toward LkHa 225 we also 
observe a dissociative shock, presumably located closer to the 
outflow source than the non-dissociative shock. In the line of 
sight towards LkHa 225 we observed absorption features due 
to solid water ice and amorphous silicates, and due to gas-phase 
H2O, CO and CO2. No solid CO2 was detected towards LkHa 
225, making this the first line of sight where the bulk of the 
CO2 is in the gas-phase. 

Key words: Shock waves - Circumstellar matter - Stars: For- 
mation - H II regions - ISM: Jets and outflows- Infrared: Stars 



1. Introduction 

The BDh-40°4124 region consists of a few tens of very 
young stars associated with the probable Herbig Ae/Be stars 
BDh-40°4124 (= VI 685 Cyg), LkHa 224 (= VI 686 Cyg) and 
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the young multiple system LkHa 225 (= V1318 Cyg), all 
of which have infrared excesses (Strom et al. 1972a; Hillen- 
brand et al. 1995). Together with NGC 6910, NGC 6914, the 
BDh-41°3731 region and IC 1318, it is part of the giant star 
forming region 2 Cyg, at a distance of about 1 kpc (Shevchenko 
et al. 1991). The stars in the BDh-40°4124 region are signifi- 
cantly younger (ages less than 1 million years) than those in 
the surrounding OB associations, with the low- and the high- 
mass stars having formed nearly simultaneously, leading some 
authors to suggest that star formation in this association might 
have been induced by the propagation of an external shock 
wave into the cloud core (Shevchenko et al. 1991; Hillenbrand 
et al. 1995). The IRAS and AFGL surveys showed a powerful 
infrared source within the BDh-40°4124 group, but the posi- 
tional uncertainty did not allow to assign the flux to one of the 
individual objects in the region. More recent submm contin- 
uum maps of the region (Aspin et al. 1994; Di Francesco et 
al. 1997; Henning et al. 1998) clearly peak on LkHa 225, sug- 
gesting that it is also dominant in the far-infrared and may be 
identified with IRAS 20187H-41 1 1 (AFGL 2557). 

Aspin et al. (1994) have shown LkHa 225 to be a triple 
system oriented north-south. The most northern and southern 
components, separated by 5", appear stellar and are photo- 
metrically variable at optical and near-infrared wavelengths. 
The middle component exhibits strong [S ll] emission, sug- 
gesting it is a nebulous knot of shock-excited material, not un- 
like many Herbig-Haro objects (Magakyan & Movseyan 1997). 
Whereas BDh-40°4124 dominates the association in the optical, 
the southern component of LkHa 225 is brightest in the mid- 
infrared (Aspin et al. 1994). A K-band spectrum of LkHa 225- 
South by the same authors shows strong ro-vibrational emis- 
sion lines of molecular hydrogen. Aspin et al. (1994) estimate 
the total luminosity of this object to be « 1600 Lq, which 
places it in the luminosity range of the Herbig Ae/Be stars. 

PaUa et al. (1994) obtained near infrared, CO, C^^O, CS, 
and H2O maser observations of the group. Their high resolu- 
tion VLA data show that a H2O maser source is clearly asso- 
ciated with LkHa 225-South. Moreover, a density concentra- 
tion in the molecular cloud (as evidenced by CS J=5^ emis- 
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Fig. 1. Combined SWS/LWS full grating spectra of (from top to bottom) BD+40°4124, LkHa 224 and LkHa 225 with the most 
prominent features identified. Apparent discontinuities in the spectra at 12.1, 27.5, 29.2 and 45.5 /xm are caused by the change to 
a bigger aperture size and may reflect the presence of extended far-infrared emission. 



sion) and a CO outflow are both associated with LkHa 225. 
In their model, LkHa 225 is at the center of a dense molecu- 
lar core of mass ss 280 M©, while BDh-40°4124 lies near the 
periphery. Continuum radio emission has been detected from 
BDh-40°4124 and a source 42" to the east of BDh-40°4124 
(Skinner et al. 1993). 



The infrared emission-line spectrum of young stellar ob- 
jects (YSOs) such as those in the BDh-40°4124 group is domi- 
nated by the interaction of the central object with the remnants 
of the cloud from which it formed. If a young stellar object pro- 
duces intense UV radiation, either due to accretion or because 
of a high temperature of the central star, an H ii region will 
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develop, producing a rich ionic emission line spectrum. At the 
edge of the H ii region, a so-called photodissociation or photon- 
dominated region (PDR) will be created, in which neutral gas 
is heated by photoelectrically ejected electrons from grain sur- 
faces. Cooling of this neutral gas occurs mainly through emis- 
sion in atomic fine-structure and molecular lines, allowing us 
to probe the physical conditions in the PDR through the study 
of its infrared emission-line spectrum. If the YSO has a strong, 
often coUimated, outflow, it will cause a shock wave as the out- 
flow penetrates in the surrounding molecular cloud, heating the 
gas. In case the shock has a velocity smaller than 40 km s~^, 
it will not dissociate the molecular material and mainly cool 
through infrared molecular transitions. Because the physical 
conditions from the pre-shock gas to the post-shock gas change 
graduaUy within such a non-dissociative shock, they are of- 
ten called C(ontinuous)-shocks. For higher shock velocities, 
the molecules will be dissociated and cooling in the shock oc- 
curs mainly through atomic and ionic emission Unes. In such a 
shock, the physical conditions from the pre- to the post-shock 
gas will change within one mean free path and therefore they 
are usuaUy termed J(ump)-shocks. In the post-shock gas of a 
J-shock, molecules will re-form and cool down the gas further. 
Therefore a dissociative shock will produce an infrared spec- 
trum consisting of a combination of ionic and molecular emis- 
sion lines. 

In this paper we present Infrared Space Observatory (ISO; 
Kessler et al. 1996) spectroscopic data obtained at the posi- 
tions of BDh-40°4124, LkHa 224 and LkHa 225. Some of the 
H2 lines included in this data-set were already analyzed in a 
previous letter (Wesselius et al. 1996). Here we will re-analyze 
an extended set of H2 data using the latest ISO calibrations, in 
combination with new data on infrared fine-structure lines, gas- 
phase molecular absorption bands and solid-state emission fea- 
tures in the BDh-40°4124 group. We will show that the infrared 
emission lines arise from the combination of a H II region and 
a PDR at the position of BD-h40°4124, from a non-dissociative 
shock at the position of LkHa 224, and from the combination 
of a non-dissociative and a dissociative shock at the position 
of LkHa 225. We will also discuss the unusual absorption-line 
spectrum seen in the line of sight towards LkHa 225 and briefly 
discuss its implications for the physical and chemical evolution 
of hot cores. 



2, Observations 

We obtained ISO Short Wavelength (2.4-45 jim) Spectrome- 
ter (SWS; de Graauw et al. 1996a) and Long Wavelength (43- 
197 jim) Spectrometer (LWS; Clegg et al. 1996) grating scans 
at die positions of BDh-40°4124, LkHa 224 and LkHa 225. For 
all three positions, both SWS scans covering the entire SWS 
wavelength range ("AOT 01") and much deeper scans covering 
small wavelength regions around particularly useful diagnostic 
Unes ("AOT 02") were obtained. For LWS only scans covering 
the entire wavelength region were obtained. A full log of the 
observations is given in Table 1. 



Table 1. Log of ISO observations in the BD-h40°4124 region. 



Object 


PsyjV 


Rev. 


Date 


JU— Z4DUUUU 


BD+40°4124 


SWS 02 


142 


07/04/1996 


181.037 




SWS 02 


159 


24/04/1996 


197.635 




SWS 01 


355 


05/11/1996 


393.283 




LWS 01 


768 


23/12/1997 


805.525 


LkHa 224 


SWS 02 


142 


07/04/1996 


180.960 




SWS 01 


858 


22/03/1998 


894.920 


LkHa 225 


SWS 02 


142 


07/04/1996 


180.996 




SWS 02 


355 


05/11/1996 


393.306 




SWS 01 


858 


22/03/1998 


894.889 




LWS 01 


142 


07/04/1996 


181.017 



Data were reduced in a standard fashion using caHbration 
files corresponding to ISO off-Une processing software (OLP) 
version 7.0, after which they were corrected for remaining 
fringing and glitches. To increase the S/N in the final spec- 
tra, the detectors were aligned and statistical outliers were re- 
moved, after which the spectra were rebinned to a lower spec- 
tral resolution. Figure 1 shows the resulting ISO spectra. Plots 
of all detected lines, rebinned to a resolution A/ A A of 2000 
with an oversampling factor of four (SWS), or averaged across 
scans (LWS), are presented in Figs. 2-A. Line fluxes for de- 
tected lines and upper limits (total flux for line with peak flux 
3 a-) for the most significant undetected Unes are Usted in Ta- 
ble 2. 

For BD-h40°4124 and LkHa 224, all detected lines appear 
at their rest wavelength within the accuracy of the ISO spec- 
trometers. For LkHa 225, all lines in the SWS spectra are sys- 
tematically shifted by « —100 km s~^. This could either be a 
reflection of the real spatial velocity of LkHa 225 or be caused 
by a slight offset between our pointing, based on the optical 
position of the source, and the position of the infrared source. 
For all three sources, the lines are unresolved, indicating that 
the velocity dispersion is smaller than a few hundred km s~^. 

For each complete spectral scan, the SWS actually makes 
twelve different grating scans, each covering a small wave- 
length region ("SWS band"), and with its own optical path. 
They are joined to form one single spectrum (Fig. 1). Because 
of the variation of the diffraction limit of the telescope with 
wavelength, different SWS bands use apertures of different 
sizes. For a source that is not point-like, one may therefore see a 
discontinuity in flux at the wavelengths where such a change in 
aperture occurs. This effect can indeed be seen in some of our 
spectra, indicating the presence of extended far-infrared emis- 
sion throughout the region. 

Since both grating spectrometers on board ISO use aper- 
tures that are fairly large compared to the separation of sources 
in most star forming regions, some caution is appropriate in 
interpreting such measurements. We created a plot with the 
positions of the SWS apertures, overlaid on a K'-band im- 
age of the region (Hillenbrand et al. 1995), shown in Fig. 5. 
As can be seen from this figure, only one strong near-infrared 
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Fig. 2. Detected lines in BD+40°4124, rebinned to heliocentric velocities. The grey Unes show the instrumental profiles for a 
point source and an extended source filling the entire aperture. 



source is included in the BD+40°4124 and LkHo 224 measure- 
ments, whereas the LkHa 225 measurement only includes the 
question-mark shaped infrared nebula formed by the point-like 
objects LkHa 225-North and LkHa 225-South and the bridg- 
ing nebula LkHa 225-Middle. ISO-LWS has a beam that is 
larger (61-83" FWHM) than the separation of the objects stud- 
ied here. Therefore the LWS measurements of BDh-40°4124 
and LkHa 225 partly cover the same region. The LWS spec- 
trum of LkHa 224 shown in Fig. 1 is a weighted average of 
these spectra. 

Accuracies of the absolute flux calibration in the SWS 
spectra range from 7% in the short-wavelength (< 4.10 i^m) 
part to « 30% in the long wavelength (> 29 /im) part (Leech 
et al. 1997). The LWS absolute flux calibration is expected to 
be accurate at the 7% level (Trams et al. 1997). Note that the 
appUed flux calibration is based on the assumption that we are 
looking at a point source. For an extended source, the diffrac- 
tion losses will be underestimated and in particular at the long- 



wavelength part of the LWS, these corrections will exceed the 
quoted uncertainty. 



3. Spectral energy distributions 

A Spectral Energy Distribution (SED) for BDh-40°4124 was 
constructed by combining UV, optical, infrared, submm and ra- 
dio data of BDh-40°4124 from literature (Wesselius et al. 1982; 
Strom et al. 1972a; Terranegra et al. 1994; Hillenbrand et al. 
1992, 1995; Cohen 1972; Lorenzetti et al. 1983; Di Francesco 
et al. 1997; Corcoran & Ray 1998; Skinner et al. 1993; Bertout 
& Thum 1982) with our new ISO spectra (Fig. 6a). As can be 
seen from this plot, the optical to submm photometry forms 
a smooth curve. At radio wavelengths the slope of the energy 
distribution appears significantly flatter. We interpret this SED 
as the sum of three components: continuum emission from the 
central star in the UV and optical, thermal emission from dust 
heated by the central star in the infrared to submm and free-free 
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Fig. 3. Same as Fig. 2 for LkHa 224. 
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emission from the stellar wind or an H ii region at cm wave- 
lengths. The SWS spectrum of BDh-40°4124 shows a rising 
continuum component starting around 13 /xm. This component 
is not present in the infrared photometry by Cohen (1972). This 
could be due to the smaller aperture size or the chopping throw 
used for background subtraction by Cohen. We interpret this 
second dust component, only visible in the ISO data, as diffuse 
emission throughout the star forming region. 

There is general consensus in the literature that the spectral 
type of BDh-40°4124 is about B2 Ve-nsh (Merrill et al. 1932; 
Herbig 1960; Strom et al. 1972b; Swings 1981; Finkenzeller 
1985; Hillenbrand et al. 1995), corresponding to an effective 
temperature Teff = 22,000 K and surface gravity logg = 4.0 
(Schmidt-Kaler 1982). With this spectral type and the optical 
photometry by Strom et al. (1972a) we arrive at a value of 
E{B — V)= (f^97 towards the optical star. This value agrees 
well with the E{B — V) derived from the diffuse interstellar 
band at 5849A (Oudmaijer et al. 1997). With this E{B - V) 
value and assuming a normal interstellar extinction law (i.e. 
Rv = Av/E{B — V)= 3.1), we corrected the data for extinc- 
tion. 

A Kurucz (1991) stellar atmosphere model with Tcff and 
log (/ corresponding to the star's spectral type (assuming solar 
abundance) was fitted to the extinction-corrected UV to optical 
SED. The fit is also shown in Fig. 6a. As can be seen from this 
figure, the SED fit is not perfect in the UV. Possibly a UV ex- 
cess due to a strong stellar wind is present in BD-h40°4124. The 
fitted Kurucz model was used to compute an apparent stellar 
luminosity L/{'iiT(f ) for BDh-40°4124, which was converted 
to an absolute stellar luminosity using the BD-i40°4124 dis- 
tance estimate of 1 kpc by Shevchenko et al. (1991). Infrared 
luminosities were computed by fitting a spline to the infrared 
to mm data and subtracting the Kurucz model from this fit. 
The derived values for the optical and infrared luminosity of 
BDh-40°4124 are 6.4 x lO^ and 3.1 x 10^ Lq, respectively. 
These luminosities place the star close to the zero-age main 
sequence position of a 10 star in the Hertzsprung-Russell 
diagram (HRD). 



We followed the same procedure to construct a SED for 
LkHa 224, combining optical, infrared and submm photometry 
from hterature (Strom et al. 1972a; Shevchenko, personal com- 
munication; Hillenbrand et al. 1995; Cohen 1972; Di Francesco 
et al. 1997) with our new ISO spectra. It is shown in Fig. 6b. 
Since LkHa 224 is strongly variable in the optical, the data 
used in the SED were selected to be obtained near maximum 
system brightness, when available. The optical to infrared pho- 
tometry of LkHa 224 forms a smooth curve in the SED, which 
we interpret as due to continuum emission from the star in the 
optical, and thermal emission from dust heated by the central 
star in the infrared. 

There is no agreement in the literature on the spectral type 
of LkHa 224. Wenzel (1980) classified the star spectroscopi- 
cally as B8e, whereas Hillenbrand et al. (1995) obtained a spec- 
tral type of B5 Ve. Corcoran & Ray (1997) derived a spectral 
type of AO for LkHa 224, based upon the relative strength of 
the He I A4471 and Mgll A4481 lines. In the optical, LkHa 
224 shows large-amplitude photometric variations (> 3™) due 
to variable amounts of circumstellar extinction (Wenzel 1980; 
Shevchenko et al. 1991, 1993). In the group of Herbig Ae/Be 
stars, these large-amplitude variations are only observed when 
the spectral type is AO or later (van den Ancker et al. 1998a). 
Therefore we favour a later spectral type than B for LkHa 224. 
The detection of the Ca ii H and K lines in absorption (Maga- 
kyan & Movseyan 1997) also seems to favour a later spectral 
type than mid-B. Recent optical spectroscopy of LkHa 224 by 
de Winter (personal communication) shows strong, variable, P- 
Cygni type emission in Ha. The higher members of the Bahner 
series (up to HI 5) are present and in absorption, together with 
the Ca II H and K lines. Numerous shell lines are present in 
the spectrum as well, which could hinder spectral classification 
solely based upon red spectra. From the available optical spec- 
troscopy, we classify LkHa 224 as A7eH-sh, corresponding to 
Teff = 7850 K. We estimate the uncertainty in this classification 
to be around 2 subclasses, or 400 K. 

With the new spectral type of A7eH-sh and the optical pho- 
tometry by Strom et al. (1972a), we derive an E{B — V) at 
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Fig. 4. Same as Fig. 2 for LkHa 225. 



maximum system brightness for LkHa 224 of (K°96. Again fit- 
ting a Kurucz model to the extinction-corrected photometry, we 
derive values for the optical and infrared luminosity of LkHa 
224 of 1.1 X 10^ and 2.7 x 10^ Lg, respectively. The good 
fit of the SED to the Kurucz model with Tes = 7850 K shows 
that the photometry is consistent with our adopted spectral type 
of A7e for LkHa 224. Comparing the derived parameters of 
LkHa 224 with the pre-main sequence evolutionary tracks by 
Palla & Stabler (1993), we see that it is located in the HRD at 
the location of a 3 x 10^ yr old 3.5 Mq star. This brings the 
age of the star in good agreement with those of the low-mass 
members of the BDh-40°4124 group (Hillenbrand et al. 1995). 

Again the same procedure as above was followed to con- 
struct a SED for LkHa 225 (Fig. 6c), combining optical to ra- 
dio measurements from literature (Shevchenko, personal com- 
munication; Hillenbrand et al. 1995; Cohen 1972; Weaver & 
Jones 1992; Di Francesco et al. 1998; Henning et al. 1998; 
Cohen et al. 1982) with our ISO data. These data refer to the 
sum of all components in this object. Since LkHa 225 shows 



strong variability in the optical and infrared, only data ob- 
tained near maximum brightness was used in the SED. Like 
in BDh-40°4124, the SED of LkHa 225 consists of a smooth 
continuum ranging from the optical to the submm. In view 
of the knowledge of the nature of LkHa 225 from literature, 
this SED must be interpreted as the sum of the three compo- 
nents of LkHa 225 in the optical to submm, each consisting 
of a hot component with circumstellar dust emission. The plot- 
ted IRAS fluxes of IRAS 20187H-4111 (Weaver & Jones 1992) 
are smaller than those in the LkHa 225 LWS spectrum. Most 
likely, this is due to the background subtraction applied on the 
IRAS data. This means that a cool component is present as well 
throughout the region, covering an area that is larger than the 
size of the IRAS detectors. In view of submm maps of the re- 
gion (Aspin et al. 1994; Henning et al. 1998), this seems plau- 
sible. 

As can be seen in Fig. 6c, below 8 /xm the flux levels in 
the SWS spectrum are considerably smaller than those from 
ground-based photometry at maximum brightness. At longer 



Table 2. Observed and extinction-corrected line fluxes and model predictions (in 10 W m ) for programme stars. 



Line 


A 

r ml 


AOT 


Beam 
[IQ-* sr] 




BD440°4124 






LkHa 224 






LkHa 225 




Observed 


Ext. corr. 


Tvlodel 


Observed 


Ext. Corr. 


Tvlodel 


Observed 


Ext. Corr. 


M!odel 


H2 0-0 S(0) 


28.2188 


S02 


1.64 


< 0.66 


< 0.68 


0.02 


< 0.95 


< 0.98 


0.03 


< 1.38 


< 1.63 


0.04 


H2 0-0S(1) 


17.0348 


S02 


1.15 


0.54±0.18 


0.57±0.19 


0.54 


1.01±0.33 


1.07±0.35 


1.04 


0.78±0.26 


1.05±0.35 


1.34 


H2 0-0 S(2) 


12.2786 


S02 


1.15 


0.71±0.23 


0.76±0.25 


0.46 


1.27±0.42 


1.36±0.45 


1.13 


2.18±0.93 


3.08±1.31 


1.79 


Hj 0-0 S(3) 


9.6649 


S02 


0.85 


0.77±0.24 


0.91±0.28 


0.90 


1.55±0.31 


1.60±0.37 


2.06 


2.80±0.73 


6.59±1.73 


5.47 


H2 0-0 S(4) 


8.0251 


S02 


0.85 


0.10±0.05 


0.11 ±0.05 


0.16 


0.69±0.23 


0.73±0.24 


0.79 


1.71 ±0.46 


2.27±0.60 


2.29 


H2 0-0 S(5) 


6.9095 


S02 


0.85 


0.25±0.11 


0.25±0.12 


0.13 


1.17±0.25 


1.21 ±0.26 


0.99 


4.03±0.63 


4.73±0.74 


4.32 


H2 0-0 S(6) 


6.1086 


S02 


0.85 


< 0.80 


< 0.83 


0.01 


< 0.63 


< 0.65 


0.08 


< 2.56 


< 3.08 


0.56 


H2 0-0 S(7) 


5.5112 


S02 


0.85 


< 0.81 


< 0.84 


0.00 


< 1.30 


< 1.35 


0.04 


3.51±0.79 


4.33±0.97 


0.45 


H2 0-0 S(8) 


5.0531 


SOI 


0.85 


< 1.83 


< 1.92 


0.00 


< 3.55 


< 3.72 


0.00 


0.81±0.40 


1.03±0.51 


0.06 


H2 0-0 S(9) 


4.6946 


S02 


0.85 


< 1.16 


< 1.22 


0.00 


< 2.55 


< 2.69 


0.00 


3.04±1.01 


3.97±1.31 


0.05 


H2 0-0S(10) 


4.4099 


S02 


0.85 


< 0.42 


< 0.44 


0.00 


< 0.42 


< 0.44 


0.00 


< 0.63 


< 0.85 


0.00 


H2 O-OS(ll) 


4.1813 


S02 


0.85 


< 0.74 


< 0.78 


0.00 


< 1.08 


< 1.15 


0.00 


0.67±0.22 


0.93±0.31 


0.00 


Ht 1-0 on"! 

jn2 1 u \iyL) 


2.4066 


SOI 


0.85 


< 0.83 


< 0.99 


0.00 


< 1.63 


^ 1 93 


0.00 


1.21±0.20 


2,98-1-0.50 


2.59 


H-1 1 -0 Cl(7'^ 
jn2 1 yj 


2.4134 


SOI 


0.85 


< 0.56 


< 0.67 


0.00 


< 1.51 


<^ 1 79 


0.00 


< 1.45 


< 3.53 


0.66 




2.4237 


SOI 


0.85 


< 0.81 


< 0.97 


0.00 


< 1.43 


< 1.70 


0.00 




3 28±0 39 


1.31 


H2 1-0 Q(4) 


2.4475 


SOI 


0.85 


< 0.52 


< 0.61 


0.00 


< 1.19 


< 1.41 


0.00 


< 1.06 


< 2.51 


0.22 


H2 1-0 0(2) 


2.6269 


SOI 


0.85 


< 2.82 


< 3.26 


0.00 


< 1.71 


< 1.98 


0.00 


< 1.07 


< 2.28 


0.70 


H2 1-0 0(3) 


2.8025 


SOI 


0.85 


< 0.63 


< 0.71 


0.00 


< 1.17 


< 1.33 


0.00 


1.75±0.15 


3.39±0.30 


2.45 


H2 1 -0 0(4) 


3.0039 


SOI 


0.85 


< 0.17 


< 0.19 


0.00 


< 0.54 


< 0.61 


0.00 


0.25±0.10 


0.44±0.18 


0.70 


H2 1-0 0(5) 


3.2350 


SOI 


0.85 


< 0.27 


< 0.30 


0.00 


< 0.89 


< 0.99 


0.00 


1.00±0.08 


1.65±0.14 


0.88 


H2 1-0 0(6) 


3.5008 


SOI 


0.85 


< 0.47 


< 0.51 


0.00 


< 0.72 


< 0.79 


0.00 


< 0.35 


< 0.53 


0.11 




185.999 


LOl 


12.2 


2 ni -1-0 32 


2 ni -1-0 32 


2.01 








ri ss-l-n 46 


ri sQ-l-n 46 


3.29 


CO J=15-14 


173.631 


LOl 


12.2 


1.95±0.29 


1.95±0.29 


1.95 


_ 


_ 


_ 


2.91±0.39 


2.92±0.39 


2.85 


CO J=16-15 


162.812 


LOl 


11.6 


1.78±0.31 


1.78±0.31 


1.85 


- 


- 


- 


2.45±0.45 


2.46±4.56 


2.41 


CO J=17-16 


153.267 


LOl 


1 1.6 


1.65±0.33 


1.65 ±0.33 


1.72 








2.01±0.35 


2.01±0.36 


1.98 


HIPfa 


7.4600 


S02 


0.85 


0.46±0.15 


0.48±0.16 


0.48 


< 0.33 


< 0.34 


- 


< 1.18 


< 1.38 


- 


HI Pf/3 


4.6539 


S02 


0.85 


0.81±0.29 


0.85±0.31 


1.21 


< 0.87 


< 0.92 




< 1.23 


< 1.61 




HI Pf7 


3.7406 


S02 


0.85 


n sQ-t-n 1 8 




1.31 


< 0.21 


< 0.22 




< 0.27 


< 0.39 




HIPf<5 


3.2971 


SOI 


0.85 


1.38±0.24 


1.51±0.25 


1.41 


< 0.58 


< 0.64 




< 0.29 


< 0.46 




HIPfe 


3.0393 


SOI 


0.85 


0.94±0.12 


1.05±0.13 


1.32 


< 0.62 


< 0.69 


- 


< 0.26 


< 0.45 


- 


HI Bra 


4.0523 


S02 


0.85 


2.71±0.54 


2.90±0.58 


2.90 


< 1.01 


< 1.08 


- 


< 1.09 


< 1.54 


- 


HI Br/3 


2.6259 


SOI 


0.85 


3.40±0.30 


3.94±0.35 


5.05 


< 1.70 


< 1.96 


- 


< 1.15 


< 2.44 


- 


rpell (^Da-^D^) 


24.0423 


S02 


1.15 


< 0.17 


< 0.18 


0.00 


< 0.15 


< 0.15 




< 0.58 


< 0.72 


0.03 


[Fell (^''_D5-^''_Dt) 


34.7133 


S02 


2.01 


< 0.40 


< 0.41 


0.00 


< 0.49 


< 0.50 




< 1.45 


< 1.63 


0.01 


fFelll C^Fq /i-'^Ft /i) 


17.9410 


S02 


1.15 


< 0.10 


< 0.10 


0.00 


< 0.19 


< 0.20 


0.00 


< 0.16 


< 0.23 


0.51 


[Fell] (^Dq/2-^Dj/2) 


25.9882 


S02 


1.15 


0.07±0.03 


0.07±0.03 


0.07 


< 0.22 


< 0.23 


0.00 


0.33±0.17 


0.40±0.20 


3.27 


[Fell] (*r'7/2-*£'5/2) 


35.3491 


S02 


2.01 


< 0.44 


< 0.45 


0.01 


< 0.77 


< 0.79 


0.00 


< 0.92 


< 1.03 


0.90 


[SI|(3p2-^Pl) 


25.2490 


S02 


1.15 


< 0.14 


< 0.14 


0.00 


0.16±0.05 


0.17±0.05 


0.16 


0.58±0.16 


0.70±0.20 


0.78 


[Sill] (-P|/2-^J'3/2) 


34.8140 


S02 


2.01 


2.32±0.77 


2.37±0.78 


2.31 


2.42±0.80 


2.48±0.82 


0.00 


3.05±1.01 


3.43±1.13 


3.06 


[01] (V2-V1) 


63.1840 


LOl 


16.3 


70.36±4.95 


70.99±4.99 


87.8 








65.97±4.78 


69.03±5.01 


58.5 


[01] (Vi-Vo) 


145.526 


LOl 


8.82 


4.15±0.40 


4.16±0.40 


3.69 








< 9.48 


< 9.57 


0.82 


[ora](3p2-3pi) 


51.8145 


LOl 


15.7 


7.99±1.20 


8.09±1.21 


8.09 








< 12.1 


< 12.9 




[ora](3pi-3po) 


88.3562 


LOl 


14.8 


8.31±0.78 


8.35±0.78 


8.35 








8.27±0.94 


8.46±0.96 




[Cn] (2pi/2-V3/2) 


157.741 


LOl 


11.6 


51.75±3.64 


51.83±3.64 


51.8 








62.54±4.63 


63.03±4.67 


0.45 



8 



M.E. van den Ancker et al.: ISO Spectroscopy of Young Intermediate-Mass Stars in the BD+40°4124 Group 



15" 


. 1 1 1 1 1 
- 


1 


1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 . 
* 


220" 










d 

45" 
w 




i 


• /' 


















1^ ' — 










If 


l\ \ *' 


41°2l'l5" 


* 

■ 


, , , 


, , , 1 


'''' 1 — 1 ■ 

0.1 pc 



20''20'"32" 31^ 30'' 29^^ 2B^ 27'' 26'' 

a (2000.) 

Fig. 5. SWS aperture positions for our measurements of (from 
right to left) BD+40°4124, LkHa 224 and LkHa 225 superim- 
posed on a K'-band image of the region. The rectangles indicate 
the apertures (in increasing size) for SWS bands 1A-2C (2.4- 
12.0 /xm), 3A-3D (12.0-27.5 /xm), 3E (27.5-29.5 /xm) and 4 
(29.5^0.5 /xm). SoUd and dashed rectangles give the orien- 
tation of the SWS aperture for the AOT S02 and SOI obser- 
vations, respectively. For clarity, only the smallest aperture is 
shown for the AOT S02 observations. 



wavelengths they agree. By extrapolating our ISO SWS spec- 
trum, we deduce that the total K band magnitude of LkHa 225 
was about 8™3 at the time of our observation. This is within 
the range of values present in the literature. Above 10 /im, the 
SWS flux levels are in agreement with ground-based measure- 
ments, suggesting that the large-ampUtude variability of LkHa 
225 is limited to the optical to near-infrared spectral range. 

The spectral type of LkHa 225 is very uncertain. Hillen- 
brand et al. (1995) give a spectroscopic classification of mid 
A-Fe for both LkHa 225-North and LkHa 225-South, whereas 
Wenzel (1980) derives a spectral type of G-Ke for the sum of 
both components. Aspin et al. (1994) suggested the southern 
component to be of intermediate-mass, based on its high bolo- 
metric luminosity. Unhke for BDh-40°4124 and LkHa 224, this 
makes the direct derivation of E{B—V) very uncertain. Adopt- 
ing E{B - V) = 4^97, corresponding to Ay = 15™4 derived 
from the silicate feature in the next section, results in an op- 
tical slope in the energy distribution that is steeper than the 
Rayleigh- Jeans tail of a black body. This could be because the 
extinction at the time of the optical observations was smaller 
than during the ISO observations. Therefore we used a rather 
uncertain value of E{B — V) = 2™36, based on the observed 
{B — V) colour of LkHa 225 and the intrinsic colours of a FO 
star, to correct the photometry of LkHa 225 for extinction. 

To obtain a rough measure for the total luminosity of LkHa 
225, we fitted a Kurucz model with Teff = 7200 K, correspond- 




BD+40"4la4 
Av = 3.01 
T,„ = 22000 K 
log g = 4.0 



LkHa 224 
A, = 2.98 
T,„ = 7850 K 
log g = 4.3 



+ 



LkHa 225 (total) 
Av = 7.32 
Ta, = 7200 K 
log g = 4.3 



Fig. 6. Observed (squares) and extinction-corrected (circles) 
spectral energy distributions of BDh-40°4124 (top), LkHa 224 
(middle) and LkHa 225 (bottom). Also shown (lines) are 
the ISO spectra and a Kurucz model fitted to the extinction- 
corrected UV to optical photometry for each object. 

ing to a spectral type of FO, to the extinction-corrected optical 
photometry. This results in values of 5.9 x 10^ and 1.6 x 10^ Lq 
for the total optical and infrared luminosities of LkHa 225. We 
conclude that at least one of the stars in the LkHa 225 system 
must be a massive (5-7 Mq) pre-main sequence object in order 
to explain these luminosities. 

4. Solid-state features 

The ISO fuU grating scans of BDh-40°4124 consist of a rel- 
atively smooth continuum, with a number of strong emission 
lines superimposed. The familiar unidentified infrared (UIR) 
emission features at 3.3, 6.2, 7.6, 7.8, 8.6 and 11.3 /xm, of- 
ten attributed to polycyclic aromatic hydrocarbons (PAHs), are 
prominently present. Possibly the 12.7 /im feature also ob- 
served in other sources showing strong PAH emission (Bein- 
tema et al. 1996) is present as well. No 9.7 /um amorphous 
silicate feature, either in emission or in absorption, is visible 
in the BDh-40°4124 SWS spectrum. The slight curvature in 
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the ground-based 8-13 /im BDh-40°4124 spectrum obtained by 
Rodgers & Wooden (1997), which they interpreted as a silicate 
feature in emission, is absent in our SWS data. 

The continuum in BDh-40°4124 seems to consist of two 
distinct components, one from 2.4-13 /xm and the other from 
13-200 /Ltm. In the SWS range, it is well fit by the superposi- 
tion of two blackbodies of 320 and 100 K. For the LWS spec- 
trum a significant excess at the long-wavelength part remains 
after the black-body fit. Probably a range in temperatures is 
present starting from 100 K down to much lower temperatures. 
The jumps in the spectrum at the positions corresponding to 
a change in aperture suggest that this component is caused by 
diffuse emission in the star forming region. 

Remarkably, the fine flux measured in the 3.3 /xm R\H fea- 
ture of BD-h40°4 124 (1.1 X lO-i^Wm-2) is much larger than 
the 2.0 X 10"'^''' W m~^ upper limit obtained by Brooke et 
al. (1993). R\H features also appear absent in the 8-13 /xm 
HIFOGS spectrum shown by Rodgers & Wooden (1997). This 
could either be an effect of the different aperture size used 
(20" X 14" for SWS versus a 1'.'4 circular aperture by Brooke 
et al.) in case the PAH emission is spatially extended, or be due 
to time-variability of the R\H emission. If the first explanation 
is correct, the PAH emission should come from a region with 
a diameter larger than 3'.'2 and hence be easily resolvable with 
ground-based imaging. 

The ISO spectrum of LkHa 224 consists of one single 
smooth continuum, with a few emission fines superimposed. 
No solid-state emission or absorption is apparent in the SWS 
spectrum. Again the sfight curvature seen in the HIFOGS 8- 
13 /Ltm spectrum of LkHa 224 (Rodgers & Wooden 1997) is 
absent in our SWS data. In view of the fact that most Herbig 
Ae stars show strong silicate emission in the 10 /xm range (e.g. 
van den Ancker 1999), the absence of the sificate feature is 
quite remarkable. It means that the dust emission must come 
from a region which is completely optically thick up to infrared 
wavelengths and yet only produce 3 magnitudes of extinction 
towards the central star in the optical. It is obvious that such 
a region cannot be spherically symmetric, but must be highly 
flattened or have a hole through which we can seen the central 
star relatively unobscured. A dusty circumstellar disk might be 
the prime candidate for such a region. The absence of silicate 
emission in LkHa 224 might then be explained by either an 
intrinsically higher optical depth in its circumstellar disk than 
those surrounding other Herbig Ae stars, a view supported by 
its relative youth, or be caused by a nearly edge-on orientation 
of the circumstellar disk. 

The infrared continuum of LkHa 224 cannot be fit by a 
single blackbody. Rather a sum of blackbodies with a range 
of temperatures starting at several hundreds Kelvins and con- 
tinuing to temperatures of a few tens of K are required to fit 
the spectrum adequately. Again the jumps in the spectrum at 
the positions corresponding to a change in aperture are present, 
suggesting that at the longer wavelengths we are looking at dif- 
fuse emission in the star forming region. However, the slope of 
the continuum in LkHa 224 is significantly different from that 
of the cool component in BDh-40°4124. This could either be 



due to a difference in the maximum temperature or to a differ- 
ence in temperature gradient of the diffuse component between 
the two positions. 

The SWS spectrum of LkHa 225 looks very different. 
Again we have a strong, smooth continuum with many emis- 
sion fines, but now also absorption lines due to gaseous molec- 
ular material as well as strong absorption features due to solid- 
state material are present. The familiar 9.7 /xm absorption fea- 
ture due to the Si-O stretching mode in amorphous silicates is 
strong. Apart from a 20% difference in the absolute flux lev- 
els, the silicate feature looks very similar to what was found 
by IRAS LRS (Olnon et al. 1986), suggesting tiiat it is con- 
stant in time. The integrated optical depth (J T{v)dv) of the 
9.7 jim feature in the ISO SWS spectrum is 196 cm~^. The 
18 /xm feature due to the bending mode of Si-O is also present 
in absorption with an integrated optical depth of 1 1 cm~^. The 
O-H bending mode of water ice at 3.0 is present in absorption 
as well (156 cm~^). 

Using intrinsic band strengths from literature (Tielens & 
AUamandola 1987; Gerakines et al. 1995), we have converted 
the integrated optical depths of the 9.7 /itm silicate and 3.0 /xm 
H2O features to column densities. The results are a column of 
1.6 X 10^* molecules for the silicates and 7.8 x 10^^ molecules 
for H2O ice. The 2: 1 ratio of these two species is within the 
range of what is found in other lines of sight (Whittet et al. 1996 
and references therein). Curiously, CO2 ice, which invariably 
accompanies H2O ice in all other lines of sight studied by ISO 
(Whittet et al. 1996; Boogert et al. 1999; Gerakines et al. 1999) 
is absent in LkHa 225. The derived upper limit on the solid 
CO2/H2O ratio is 0.04, much less than the canonical value of 
0.15. 

No aperture jumps are visible in the SWS spectrum of 
LkHa 225. The bulk of the continuum flux must therefore 
come from an area that is small compared to the smallest SWS 
aperture (20" x 14"). The discontinuity between the SWS and 
LWS spectra shows that although LkHa 225 itself is proba- 
bly the strongest far-infrared source in the region, the total far- 
infrared flux of the group is dominated by the diffuse compo- 
nent. Like in LkHa 224, the SWS spectrum of LkHa 225 can- 
not be fit well with a single blackbody. Again a sum of black- 
bodies with a range of temperatures starting at several hundreds 
and continuing to temperatures of a few tens of K are required 
to fit the spectrum adequately. 

PAHs appear absent in our LkHa 225 SWS spectrum. In 
view of the claim by Deutsch et al. (1994) that the southern 
component of LkHa 225 is the strongest UIR emitter in the re- 
gion, this is remarkable. One possible explanation for this ap- 
parent discrepancy could be the large degree of variabiUty (> 
3™ in K\ Allen 1973) of LkHa 225. In this case the compo- 
nent responsible for the PAH emission should have been much 
fainter than the other component at the time of observation. 
However, Aspin et al. (1994) showed the southern component 
to be dominant longward of 3 jim. Comparison of our LkHa 
225 spectrum with infrared photometry by previous authors 
shows that the source was not near minimum brightness at the 
time of the SWS observation, making the explanation that the 
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Fig. 7. Detected gas-phase molecular absorption bands in LkHa 225 with the continuum normaUzed to unity (top curves). Also 
shown (bottom curves) are the computed absorption spectra, shifted for clarity. The large plots show the AOT 01 SWS data, 
whereas the insets show the AOT 02 data containing molecular bands, T and b were required to be identical in all fits, (a) 
Observed CO ro-vibrational fines, compared to a Tex = 300 K, 6 = 5 km s"\ A^(CO) = 10^° cm~2 model, (b) H2O 1/2 band, 
with a Tex = 300 K, & = 5 km s-\ A^(H20) = lO^'' cm-^ model, (c) gas-phase CO2 V3 band, with a Tex = 300 K, 6 = 5 km s-\ 
NiC02) = 3 X 10^'^ cm-2 model, (d) gas-phase CO2 1^2 band, with a Tex = 300 K, 6 = 5 km s-\ Ar(C02) = 3 x 10^'^ cm'^ 
model. 



southern component of LkHa 225 was faint improbable. The 
most likely explanation for the apparent discrepancy between 
our non-detection and the reported UlR emission (Deutsch et 
al, 1994) might therefore be that the UIR emission within the 
southern component of LkHa 225 is variable in time. Note that 
because the gas cooling timescale is much longer than the PAH 
cooling timescale, if PAH emission has been present within re- 
cent years, gas emission from a remnant PDR might be ex- 
pected in the vicinity of LkHa 225, 

Since extinction in the continuum surrounding the 9,7 
feature is small compared to the extinction within this fea- 
ture, the extinction A\ at wavelength A across a non-saturated 
9,7 fim feature can simply be obtained from the relation Ax = 
— 2.51og(///o), Using an average interstellar extinction law 
which includes the sificate feature (Fluks et al, 1994), we can 
then convert these values of Ax to a visual extinction, resulting 
in a value of Ay = 15^4 ± 0',"2 toward LkHa 225, This value 
is significantly smaller than the Ay « 50™ obtained from the 
C^^O column density towards LkHa 225 (Paha et al, 1995), 
Although smaller, our estimate of Ay is still compatible with 



the value of 25™ ± 9™ derived from the ratio of the 1-0 S(l) 
and Q(3) lines of H2 (Aspin et al, 1994), 

5. Gas-phase molecular absorption lines 

In the SWS spectrum of LkHa 225 several absorption bands 
are present due to a large column of gas-phase CO, CO2 and 
H2O in the line of sight. They are shown in Fig, 7, Molecu- 
lar absorption bands are absent in the BD+40°4124 and LkHa 
224 spectra. We compared the observed gas-phase absorption 
lines to synthetic gas-phase absorption bands computed using 
molecular constants from the HITRAN 96 database (Rothmann 
et al, 1996), First the relative population of the different levels 
within one molecule were computed for an excitation temper- 
ature Tex, For the relative abundances of the isotopes, values 
from the HITRAN database were adopted, A Voigt profile was 
taken for the lines, with a Doppler parameter b = ^^^^ ^ for 
the Gaussian component of the Voigt profile and the lifetime 
of the considered level times b/Air for the Lorentz component. 
The optical depth tx as a function of wavelength was then com- 
puted by evaluating the sum of all integrated absorption coef- 
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ficients, distributed on the Voigt profile. Synthetic absorption 
spectra were obtained by computing exp(— ta). This proce- 
dure for computing absorption spectra is nearly identical to that 
followed by Helmich et al. (1996) and Dartois et al. (1998). 
For comparison with the observations, the spectra were con- 
volved with the ISO SWS instrumental profile. For the AOT 02 
observations, the instrumental profile is nearly Gaussian with 
a wavelength-dependent resolving power R = FWHM/X 
ranging between 1400 and 3000 (Leech et al. 1997). For our 
AOT 01 observations, the resolution is a factor four lower. 

Best fits of the synthetic absorption spectra to the observed 
CO, CO2 and H2O bands are also shown in Fig. 7. Because of 
the relatively low signal to noise in the spectra, the fits for the 
different absorption bands were not done independently, but 
Tex and b were required to be the same for all species. This 
procedure is only valid if the different molecular species are co- 
spatial. The results of similar fitting procedures in other lines 
of sight (e.g. van Dishoeck et al. 1996; Dartois et al. 1998; 
Boogert 1999) suggest that this is the case. Satisfactory fits to 
the data could be obtained for excitation temperatures of a few 
100 K and a Doppler parameter b between 3 and 10 km s~^ 
The derived column of warm CO2 is several times 10^^ cm^^. 
For H2O it is about 10^^ cm~^, whereas the column of gas- 
phase CO is in the range 10^^-10^^ cm~^. 

The derived values of b are much higher than 
those expected from purely thermal broadening 
(0.1290 i/TexM km s~\ with A the molecular weight 
in amu; Spitzer 1978). This must be caused by macroscopic 
motion within the column of warm gas in our line of sight 
toward LkHa 225. The fact that our deduced values of b 
are within the range of those found toward hot cloud cores 
(Helmich et al. 1996; van Dishoeck & Helmich 1996; Dartois 
et al. 1998) suggests that it is due to turbulence within the 
cloud rather than a dispersion in large-scale motion such as 
infall or rotation. 

6. Hydrogen recombination lines 

The ISO SWS spectrum of BD-h40°4124 contains several emis- 
sion lines of the Pfund and Brackett series of H I. No H I lines 
were detected in LkHa 224 and LkHa 225. The lines fluxes for 
BD-h40°4124 were corrected for extinction using the value of 
Ay = 3™0 derived in section 3. We combined these data with 
ground-based measurements of emission Unes in the Brackett 
and Paschen series (Harvey 1984; Hamann & Persson 1992; 
Rodgers & Wooden 1992) for BDh-40°4124 to create plots of 
the H I line flux ratio as a function of upper level quantum num- 
ber n (Fig. 8). The fact that the values for the Brackett series lie 
on a continuous curve shows that they are not influenced much 
by the different beam sizes used, indicating that they must orig- 
inate in a compact region. 

In Fig. 8 we also show predicted values for Menzel & Baker 
(1938) Case B recombination for a range of temperatures from 
1000 to 30,000 K and densities from 10^ to 10^ cm'^, taken 
from Storey & Hummer (1995). The fit of the Case B values 
to the data is poor, indicating that at least the lower lines in 
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Fig. 8. Extinction-corrected H I fine flux ratios in BDh-40°4124 
for (from top to bottom) the Pfund, Brackett and Paschen series. 
Also shown (dotted lines) are theoretical values for Case B re- 
combination for T = 1000-30,000 K and n = lO^-lO^ cm'^. 
The dashed Unes give model predictions for an optically thick 
static slab of gas with temperatures (from bottom to top) be- 
tween 1000 and 100,000 K. The gray soHd line shows the best 
fit of the LTE fully ionized wind model. 



the Pfund and Brackett series are not opticaUy thin, as assumed 
in Case B recombination theory, but have r ^ 1. An error in 
our estimate of the extinction towards the line emitting region 
cannot reconcile our observations with the Case B values. The 
poor fit of the Case B model predictions to the observed H i line 
fluxes means that the lines do not originate in the low-density 
H II region responsible for the [O III] emission (section 8). 

In the limit that all lines are optically thick, the line flux ra- 
tios of a static slab of gas in local thermodynamic equilibrium 
(LTE) will be dominated by the wings of the line profile. Such a 
slab of gas will produce line fluxes Ix oc Bx (Tg) ffj^ , with Bx 
the Planck function for an electron temperature Te, and fij the 
oscillator strength for the transition (Drake & Ulrich 1980). In 
Fig. 8 we also show these values for a range of electron temper- 
atures between 1000 and 100,000 K. Again this simple model 
cannot reproduce the observed line flux ratios in BDh-40°4124, 
especially for Pfa/Bra and Br7/Bra. The most Ukely explana- 
tion for the poor fit is that the lines are more broadened than due 
to the Stark wings of the line profile, probably due to macro- 
scopic motion of the gas. 
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The poor fit of both the Case B values and the optically 
thick static slab model shows that we are not looking at an op- 
tically thick static slab of gas, but at an optically thick ionized 
gas with a significant velocity dispersion. Most likely it origi- 
nates either in an ionized wind, or in an ionized gaseous disk 
surrounding BD-i40°4124. Without information on the kine- 
matics of the gas, the distinction between these two cases is 
hard to make. To qualitatively show that such a model can in- 
deed reproduce the observations, we also show in Fig. 8 the best 
fit of a simple LTE model for a totally ionized wind following 
the treatment of Nisini et al. (1995). This best fit model has a 
rather high mass loss rate (~ 10~^ M© yr~^), arising from a 
compact region (« 12 R^) in radius. However, the assumption 
of LTE will cause this simple model to significantly overesti- 
mate the mass loss rate (Bouret & Catala 1998), whereas the 
assumption that the wind is completely ionized will cause an 
underestimate of the size of the line emitting region. There- 
fore the values given here should be approached as lower, re- 
spectively upper limits. Also, a model for an ionized gaseous 
disk surrounding BDh-40°4124 may give equally good fits to 
the data. 

If the lower lines in the Pfund and Brackett series are op- 
tically thick, the region in which these lines originate should 
also produce optically thick free-free emission at radio wave- 
lengths. In the simple wind model employed above, the Bra 
line flux is related to the 6 cm continuum flux by the relation 

(mJy) = 1.4 X 10^^ Fb™ (W m'^) v^'^'^, with V2 the wind 
velocity in units of 100 km s"'^ (Simon et al. 1983). Assum- 
ing V2 = 1, fliis relation yields 0.39 mJy for BDh-40°4124, in 
excellent agreement with the observed 6 cm flux of 0.38 mJy 
(Skinner et al. 1993). Although highly simplified, this analysis 
does show that the region responsible for the H l line emission 
can also produce sufficient amounts of radio continuum emis- 
sion to account for the observed fluxes and that our interpreta- 
tion of the radio component in the SED as free-free emission 
(section 3) is correct. 

7. Molecular hydrogen emission 

At all three observed positions in the BDh-40°4124 region we 
detected emission from pure rotational lines of molecular hy- 
drogen. They are listed in Table 1 and are shown in Figs. 2-A. 
For BDh-40°4124 and LkHa 224, 0-0 S(0) to S(5) were de- 
tected. For LkHa 225 pure rotational lines up to S(10) as well 
as several ro- vibrational lines were detected. It is clear that in 
LkHa 225 the 0-0 S(6) line escaped detection because of its 
unfortunate location in the H2O gas phase 1^2 absorption band 
(Fig. 7b). The line flux of the 0-0 8(5) line may also be sUghtly 
(sd 10%) underestimated because of this. We have detected sig- 
nificantly more lines than the preliminary results reported in the 
A&A ISO First Results issue (Wesselius et al. 1996), bofli be- 
cause of improvements in the SWS data reduction since 1996 
and because additional data were obtained. Therefore a new 
analysis of these H2 lines is in place. 

From the H2 line fluxes I{J) it is possible to calculate 
the apparent column densities of molecular hydrogen in the 
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Fig. 9. H2 excitation diagrams for BDh-40°4124 (top), LkHa 
224 (middle) and LkHa 225 (bottom). ISO observations of 
pure rotational lines are indicated by the filled dots. Squares, 
triangles and stars indicate measurements of ro-vibrational 
lines from the SWS spectra and from Aspin et al. (1994). For- 
mal errors for most measurements are smaller than the plot 
symbol. The dashed lines gives the Boltzmann disttibution fits 
to the low-lying pure rotational lines. 



upper J levels, averaged over the SWS beam, N{J), using 
N{J) = with A the wavelengfli, h Planck's con- 

stant and c the speed of Ught. The transition probabilities Aij 
were taken from Turner et al. (1977). Line fluxes were cor- 
rected for extinction using /(J) = JobslO"^'-'''-'/^'^ using the 
average interstellar extinction law by Fluks et al. (1994). For 
both BDh-40°4124 and LkHa 224 we used Ay = 3^0, valid 
for the optical stars. For LkHa 225 we used the Ay = 15™4 
derived in section 4. 

A useful representation of the H2 data is to plot the log of 
N{J)/g, the apparent column density in a given energy level 
divided by its statistical weight, versus the energy of the upper 
level. The statistical weight g is the combination of the rota- 
tional and nuclear spin components. It is (2 J + 1) for para- 
H2 (odd J) and (2 J -|- 1) times the Ortho/Para ratio for ortho- 
H2 (even J). We have assumed the high temperature equilib- 
rium relative abundances of 3: 1 for the ortho and para forms of 
H2 (Burton et al. 1992). The resulting excitation diagrams are 
shown in Fig. 9. For LkHa 225 we also included the measure- 
ments of the H2 ro-vibrational lines by Aspin et al. (1994). For 
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Fig. 10. Continuum-subtracted 2.35-3.30 SWS spectrum 
of LkHa 225. Also shown (bottom curve) is a synthetic H2 
spectrum with Texc = 710 K, corrected for an extinction of 
Av = 15^4. 



the lines in common, their values differ somewhat with ours, al- 
though the line ratios agree. This is probably an effect of their 
smaller beam size and indicates that in our measurements the 
beam filling factor is smaller than one. In Fig. 9 the values of 
Aspin et al. were scaled to our measurements. 

In case the population of the energy levels is thermal, all 
apparent column densities should lie on a nearly straight line 
in the excitation diagram. The slope of this line is inversely 
proportional to the excitation temperature, while the intercept 
is a measure of the total column density of warm gas. As can be 
seen from Fig. 9, this is indeed the case in all three stars for the 
pure rotational lines with an upper level energy up to 5000 K. 
The fact that the points for ortho and para H2 lie on the same 
line proves that our assumption on their relative abundances is 
correct. We have used the formula for the H2 column density 
for a Boltzmann distribution given by Parmar et al. (1991) to 
fit our data points in the low-lying pure rotational levels, using 
the rotational constants given by Dabrowski (1984), varying 
the excitation temperature and column density. The results are 
shown as the dashed lines in Fig. 9. 

Van den Ancker et al. (1998b and in preparation) have em- 
ployed predictions of H2 emission from PDR, J-shock and C- 
shock models (Burton et al. 1992; Hollenbach & McKee 1989; 
Kaufman & Neufeld 1996), to determine the excitation temper- 
ature from the low-lying pure rotational levels as a function of 
density n and either incident FUV flux G (in units of the aver- 
age interstellar FUV field Go = 1.2 x lO"'^ erg cm^^ s~^ sr~^; 
Habing 1968) or shock velocity Vs in an identical way as was 
done for the observations presented here. They arrived at the 
conclusion that the PDR and J-shock models allow a fairly 
small (200-540 K) range of excitation temperatures, whereas 
for C-shocks this range is much larger (100-1500 K). In the 
model predictions for shocks, Texc does not depend much on 
density, whereas for PDRs it does not depend much on G. Once 
the mechanism of the H2 emission is established, it can there- 
fore be used to constrain Vg or n in a straightforward way. 



The values for T^xc measured for LkHa 224 and LkHa 225 
fall outside the range of Texc values predicted by either PDR or 
J-shock models. They are similar to the values found in other 
regions containing shocks (van den Ancker et al. 1998b). If the 
observed pure rotational H2 emission is indeed due to a planar 
C- shock, its shock velocity should be around 20 km s~^. In 
section 4 we saw that our BDh-40°4124 SWS full grating scan 
contains strong emission in the PAH bands. This indicates that 
a PDR is present at that location. If we assume that the observed 
H2 emission is entirely due to this PDR, the observed Texc of 
470 K points to a rather high (^ 5 x 10^) density within the 
PDR. 

For LkHa 225, ro- vibrational as well as pure rotational 
emission from H2 was detected (Fig. 10). Both the higher pure 
rotational lines and the ro-vibrational Unes significantly devi- 
ate from the straight line defined by the lower pure rotational 
lines in the excitation diagram and do not form a smooth fine by 
themselves, as expected in the case of multiple thermal compo- 
nents. Such a situation can either occur because of fluorescence 
of H2 in a strong UV field (Black & van Dishoeck 1987; Draine 
& Bertoldi 1996), be due to X-ray or EUV heating (Wolfire & 
Konigl 1991; Tine et al. 1997), or be due to the formation en- 
ergy of re-formed H2 in post-shock gas (Shull & Hollenbach 
1978). Although the H2 spectra resulting from the possible for- 
mation mechanisms do differ in the relative prominence of tran- 
sitions with Av > 2 (Black & van Dishoeck 1987), the distinc- 
tion between these situations is rather difficult to make with the 
detected lines in LkHa 225. 

By summing all the apparent column densities for the ro- 
vibrational and higher pure rotational lines we obtain an esti- 
mate for the observed total column of non-thermal H2 in LkHa 
225 of 2 X 10^^ cm-2, corresponding to 2 x IQ-*^ Mq or 0.3 
earth masses within the SWS beam. This value of A^(H2) is 
within the range predicted by UV fluorescence models. How- 
ever, the absence of PAH emission in LkHa 225 at the time of 
the ISO observations suggests that there was no strong UV ra- 
diation field at that time. Therefore we consider the possibiUty 
that the ro-vibrational H2 emission originates in a shock more 
likely. 

Interestingly, in LkHa 225 we have detected two pairs of 
lines arising from the same energy level: the 1-0 Q(l) & 1- 
0(3) and flie 1-0 Q(3) & 1-0 0(5) lines. For these pairs of 
lines N{J) must be identical. Therefore they can be used to 
derive a differential extinction between the two wavelengths 

of the lines using A{\i) - A{X2) = 2.5 log (ifif^^) ■ Us- 
ing the interstellar extinction curve by Fluks et al. (1994), we 
can then convert the differential extinction to a value of Ay. 
The results of this procedure are values for Ay of 60™ ± 20™ 
and 40™ ± 20™ for tiie 1-0 Q(l) & 1-0 0(3) and 1-0 Q(3) & 
1-0 0(5) line pairs, respectively. These extinction values are 
somewhat larger than the 15™4 ± 0™2 derived from the sili- 
cate feature in section 4. But they agree well with Ay ~ 50™ 
obtained from the C^^O column density towards LkHa 225 
(Palla et al. 1995). If the slope of the infrared extinction law 
is approximately interstellar, the extinction towards the region 
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Fig. 11. CO excitation diagrams for BD+40°4124 (top) and 
LkHa 225 (bottom). The dashed lines give the fit of the Boltz- 
mann distribution to the data points. 

producing the low-lying pure rotational H2 lines in LkHa 225 
must be smaller than this value though; such a large extinction 
would cause the extinction-corrected 0-0 S(3) line to deviate 
significantly from the thermal distribution line in the excitation 
diagram, since it is located within the 9.7 /xm silicate feature 
in the extinction curve. We conclude that in LkHa 225 either 
there must be two distinct sources of H2 emission, a 710 K 
thermal component suffering little extinction - presumably a 
C-shock - and a heavily extincted, deeply embedded compo- 
nent - possibly a J-shock, or the slope of the infrared extinction 
law towards this region must be much steeper than interstellar. 

8. Carbon-monoxide emission lines 

At both the positions of BD-h40°4124 and LkHa 225, several 
emission lines due to the v = O-O transitions of gas-phase CO 
were detected in the long-wavelength part of the LWS spec- 
tra (Table 2). Similar to what was done for the H2 emission in 
the previous section, we constructed CO excitation diagrams, 
using molecular data from Kirby-Docken & Liu (1978). They 
are shown in Fig. 11. The temperature and colunm of CO re- 
sulting from the Boltzmann fit to these excitation diagrams are 
500 K and 8.1 x 10^^ cm~^ and 300 K and 1.9 x lO" cm'^ 
for BD-h40°4124 and LkHa 225, respectively. We estimate the 
errors in these fit parameters to be around 50 K in temperature 
and 50% in CO column. 

The observed CO lines have critical densities of around 
10^ cm~^. Therefore the observation of these lines also implies 
densities of at least this order of magnitude in the originating 
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Fig. 12. Behaviour of the [Oiii] 51.8 /im/88.4 /xm line flux 
ratio as a function of electron density rig (solid Unes). The 
curves for temperatures between 5,000 and 50,000 K (from 
right to left) are nearly identical. The hatched region between 
the dashed Unes shows the error interval of the observed [O iii] 
ratio in BD-h40°4124. 



region. The CO excitation temperature of 500 K found at the 
position of BDh-40°4124 is comparable to that found from the 
H2 lines. The CO/H2 mass fraction of 0.01 is compatible for 
what one would expect for the warm gas in a PDR. However, 
assuming that densities of around 10® cm""^ would exist in the 
entire BDh-40°4124 PDR would be implausible. We therefore 
conclude that if these CO lines and the H2 both originate in the 
large-scale environment of BDh-40°4124, the PDR must have a 
clumpy structure (e.g. Burton et al. 1990). 

The CO emission arising from the neighbourhood of LkHa 
225 is remarkably different from that seen in H2: here the ex- 
citation temperature of the CO is much lower than that found 
in H2. The CO/H2 mass fraction of 0.03 is also much higher 
than that found at the position of BDh-40°4124. It seems likely 
that another component than the C-shock responsible for the 
H2 emission needs to be invoked to explain the low CO ex- 
citation temperature. Again this region must have densities of 
around 10® cm^''. A possible disk or extended envelope around 
LkHa 225 could have both the temperature and densities re- 
quired to explain the observed CO spectrum. 

9. Atomic fine structure lines 

In order to interpret our observations as arising in either a PDR, 
C-Shock or J-shock, we compare our results with theoretical 
models of such regions (Tielens & HoUenbach 1985; HoUen- 
bach & McKee 1989; Kaufman & Neufeld 1996). Important 
constraints come from the observed fine structure lines. The 
mere presence of certain lines helps us to identify the mecha- 
nism responsible for the observed emission. In contrast to both 
C- and J-shocks, PDRs do not produce significant quantities 
of [S I] emission (Tielens & HoUenbach 1985). C-shocks only 
contain trace fractions of ions and hence cannot explain the 
[Fe II], [Si II] or [C ll] emission. 
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n [cm ^] 


G[Go] 


n [10-* sr] 


n [cm *] 


Vs [km s-^] 


n [10-* sr] 


n [cm ^] 


Va [km s-^] 


n [10-* sr] 


BD+40°4124 


10^/10^ 


lOVlO^ 


0.55/0.05 














LkHa 224 








« 10^ 


20 


«0.10 








LkHa 225 








« 10^ 


20 


«0.20 


10^-10^ 


60 


0.12 



In BD+40°4124 the H2 Unes and the presence of PAH 
emission strongly suggest the presence of a PDR. No [S i] was 
detected, so no shock is required to explain the emission from 
this region. The [Fe ll] and [Si ll] emission observed at the po- 
sition of BDh-40°4 124 may also originate in a PDR. The [O III] 
emission Unes detected by SWS cannot be produced in a PDR. 
The ratio of the [O III] 52 and 88 /^m lines is sensitive to the 
electron density rie in the fine forming region. It hardly de- 
pends on the temperature. We computed the expected ratio for 
a range of temperatures and densities, in the approximation that 
we can treat [O iii] as a five level ion following the treatment by 
AUer (1984), using collisional strengths from Aggarwal et al. 
(1982). The results are shown in Fig. 12. We conclude that the 
[O III] lines observed by LWS are formed in a H 11 region with 
rie = 0.6-2.2 X 10^ cm'^. In view of the 22,000 K effective 
temperature of BD-h40°4124, it is surprising that [O III] emis- 
sion is detected, since simple photo-ionization models predict 
that for stars with T^s < 32,000 K most oxygen should remain 
neutral. Clearly another ionization model must be present here. 
Possibly the region also responsible for the H i lines could col- 
hsionally ionize the oxygen and we are observing the outer, 
low-density part of an ionized wind. 

Apart from the PDR, a Hll region surrounding 
BDh-40°4124 may also contribute significantly to the ob- 
served [Fell] and [Sill] emission. Therefore we have used 
the photo-ionization code CLOUDY (version 90.04; Ferland 
1996) to generate model predictions for line strengths for an 
Hll region surrounding BDh-40°4124. A Kurucz (1991) model 
with parameters appropriate for BDh-40°4124 was taken as 
the input spectrum, and a spherical geometry and a constant 
electron density of 2 x 10^ cm-^ throughout the H II region 
were assumed. According to this model, 30% of the observed 
[Si II] 34.82 /xm flux and 40% of flie 25.99 /xm [Fe 11] flux will 
originate in the H II region. From comparing the remainder of 
the fluxes of the atomic fine structure lines to the PDR models 
of Tielens & HoUenbach (1985), we derive that the PDR 
should be fairly dense (« 10* cm-'^) and have an incident 
far-UV field of about 10"^ Gq. The extent of this PDR should 
be around 180 square arcseconds, or about half of the SWS 
beam. 

In view of the observed H2 spectrum (section 7), a C-shock 
appears the most likely candidate for explaining the observed 
[S I] emission in LkHa 224. A C-shock with a shock velocity of 
around 20 km s"^ can explain both the H2 and [S l] emission. 
However, the density and extent of the region are poorly con- 



strained: the emission may either arise in a dense (~ 10^ cm"^) 
region of only a few square arcseconds, or arise in a region of 
about 10^-10* cm"^ with a high beam filling factor. Since it 
is observed in the largest SWS aperture, the observed [Sill] 
emission could be due to contamination from LkHa 225 (see 
Fig. 5). 

For LkHa 225 the situation appears more complicated. The 
detection of [S i] clearly shows that a shock must be present 
in the region. A C-shock can explain the observed thermal H2 
component, with similar parameters as for LkHa 224, but can 
only explain part of the [S i] emission. The [Fe 11] and [Si 11] 
emission could either come from a deeply embedded H 11 re- 
gion, a deeply embedded PDR, or a J-shock. In view of the ne- 
cessity to have an additional source of [S i] emission, the most 
likely candidate might be a J-shock. A J-shock model with Vg 
w 60 km s-^ and n w 10'*-10'' cm-'^ can indeed reproduce 
the observed [S l] and [Si 11] intensities. However, to repro- 
duce the observed [Fe 11] 26.0 iim strength, iron would have 
to be about eight times more depleted than assumed in the Tie- 
lens & HoUenbach PDR models. The observed [C 11] 157.7 /xm 
emission at the position of LkHa 225 cannot be explained by 
shocks. This could very well be due to contamination by the 
BDh-40°4124 PDR or flie galactic [C 11] background, however. 

Fit parameters for the PDR, J-shock and C-shock models 
used to explain the H2, CO and fine structure emission lines 
are Usted in Table 3. The sunmied model predictions for the 
Une fluxes are listed in Table 2. 

10. Discussion and conclusions 

In the previous sections we have seen that in a small OB asso- 
ciation like the BDh-40°4124 group, different phenomena con- 
tribute to at first glance similar line spectra. The emission at 
the position of BD-h40°4124 is well reproduced by the com- 
bination of a H II region around BDh-40°4124 with a PDR. In 
Section 7 we derived values of the far-UV radiation field for 
its PDR from the observed emission lines. It is useful to com- 
pare those values to those expected from the central star to see 
whether that is a sufficient source of radiation, or another excit- 
ing source needs to be invoked. From the Kurucz (1991) model 
with Toff = 22,000 K and log 5 = 4.0, fitted to the extinction- 
corrected UV to optical SED of BDh-40°4124 (section 3), we 
compute that BD-h40°4124 has a total far-UV (6-13.6 eV) lu- 
minosity of 3.2 X 10^ L0. At a location 5 x 10^ AU from the 
central star, this field will have diluted to the value of 10^ Go 
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Table 4. Compaiison of gas-phase and solid-state absorption column densities (in cm~^) towards YSOs. 



Object 




gas-phase column 




solid colurrm 




gas/solid ratio 


Ref. 


CO 


CO2 


H2O 


CO 


CO2 


H2O 


CO 


CO2 


H2O 


LkHa 225 




1-5 X 10^'^ 


0.5-3 X 10" 


< 6 X 10" 


< 3 X 10^® 


7.8 X 10" 


> 16 


>3 


« 13 


(1) 


TTau 




< 1 X IQi'^ 


< 1 X 10^^ 


< 5 X 10" 


1.4 X 10" 


5.4 X 10" 


> 20 


< 0.7 


< 1.9 


(2) 


AFGL4176 


2 X 10^® 


5 X 10^° 


2 X 10^® 


< 5 X 10^** 


1.2 X 10" 


9 X 10" 


>400 


«0.04 


!«2.2 


(3,4,5,6) 


W33A 


4.0 X 10" 




1.0 X 10^* 


9.0 X 10" 


1.5 X 10^** 


2.8 X 10" 


w40 




wO.04 


(3,7,8,9) 


AFGL 2136 


1.0 X 10" 


1 X 10^^ 


2 X 10^* 


1.8 X 10" 


6.1 X 10" 


5.0 X 10^® 


w50 


wO.02 


wO.4 


(3,4,5,6) 


RAFGL 7009S 


6.1 X 10^* 


1.0 X 10" 


> 2 X 10^® 


1.8 X 10^* 


2.5 X 10^* 


1.1 X 10" 


«3.4 


«0.04 


>0.18 


(10) 


Elias 29 


> 1 X 10^'' 




> 9 X 10^^ 


1.7 X 10" 


6.5 X 10" 


3.0 X 10^** 


> 55 




> 0.3 


(9) 


AFGL 2591 


1.1 X 10" 


1 X 10^'' 


2 X 10^* 


< 4 X lO^'' 


2.7 X 10" 


1.7 X 10^* 


> 270 


«0.04 


« 1.1 


(3,4,5,6) 


NGC 7538 IRS9 


1.4 X 10" 


8.0 X 10^^ 


< 3 X 10" 


1.0 X 10^* 


1.2 X 10^* 
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(1996b); (5) van Dishoeck & Helmich (1996); (6) Chiar et al. (1998); (7) Gerakines et al. (1999); (8) Boogert (1999); (9) Dartois et al. (1998). 



obtained from the infrared lines. At the BDh-40°4124 distance 
of 1 kpc, this corresponds to an angular separation of 5", con- 
sistent with the constraints imposed by the SWS aperture. We 
conclude that BDh-40°4124 can indeed be responsible for the 
incident UV flux upon the PDR observed in that vicinity. Com- 
pared to the far-UV luminosity of BDh-40°4124, the other two 
sources only show little UV emission (3 and 6 Lq for LkHa 

224 and 225, respectively), explaining why we don't observe a 
PDR at those positions. 

For LkHa 224 we showed that the infrared emission- 
line spectrum could be explained by a single non-dissociative 
shock. A C-shock with very similar parameters may also be 
present at the position of LkHa 225. Comparing our ISO SWS 
aperture positions (Fig. 5) with the image of the bipolar out- 
flow centered on LkHa 225 (Palla et al. 1995), we note that 
this outflow covers both LkHa 224 and LkHa 225. One possi- 
ble identification for the C-shock in both LkHa 224 and LkHa 

225 might therefore be the shock created as the outflow orig- 
inating from LkHa 225 drives into the surrounding molecular 
cloud. 

Apart from the large-scale C-shock, we have seen that prob- 
ably a J-shock is present in LkHa 225 as well. This can be nat- 
uraUy identified with the Herbig-Haro knot LkHa 225-M ob- 
served in [S II] in the optical (Magakyan & Movseyan 1997). It 
remains unclear whether the ro-vibrational H2 lines observed 
in this object are due to H2 fluorescence by Lya photons or to 
H2 re-formation in the post-shock gas. Spectroscopy of lines 
with higher v is required to clarify this. 

In the line of sight toward LkHa 225, a large column con- 
sisting of warm gas-phase CO, CO2 and H2O and soUd water 
ice and silicates was detected. The column of gas-phase water 
was found to be about 10 times higher that that of water ice. No 
CO2 ice was detected toward LkHa 225, demonstrating that the 
bulk of the CO2 is in the gas-phase. In Table 4 we compare the 
derived gas/sohd ratio for LkHa 225 with those found in the 
lines of sight towards other young stellar objects. Both the H2O 
and CO2 gas/solid ratios are much higher than those found in 
other YSOs. For CO2, the lower limit to the gas/solid is even a 



factor 100 higher than that found in the source with the previ- 
ously reported highest ratio! At this point we can only specu- 
late as to why this is the case. One possible explanation could 
be that the radiation field of BDh-40°4124 is heating the outer 
layers of the LkHa 225 envelope, evaporating the CO2 ice. If 
this is the case, BDh-40°4124 should be located closer to the ob- 
server than LkHa 225, suggesting that LkHa 225 is located at 
the far side of the "blister" in which the BD-h40°4124 group is 
located. Another explanation could be that the LkHa 225 core 
itself is unusually warm. Since the abundance of gas-phase wa- 
ter is sensitive to the temperature (Chamley 1997), this would 
also explain the high abundance of water vapour in the line of 
sight toward LkHa 225. Because LkHa 225 is less massive 
than the luminous objects in which the previous gas/soUd ra- 
tios have been determined, this high temperature might reflect 
a more evolved nature of the LkHa 225 core. Future research, 
such as a temperature determination from gas-phase CO mea- 
surements, should be able to distinguish between these possi- 
bilities and show whether we have discovered the chemically 
most evolved hot core known to date, or are looking through a 
Une of sight particularly contaminated by a nearby OB star. 
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